Luminosity vs. Line-Width for Spiral Galaxies 


Traditionally, the most accepted measurement of the Hubble 
Constant comes from the luminosity-line width relation for spiral 
galaxies, otherwise known as the Tully-Fisher relation. Since the 
technique uses star-forming galaxies, it can be well-calibrated 
with Cepheids, and, since it involves the luminosity of entire 
galaxies, it can easily be used to measure the undisturbed Hub- 
ble Flow. Unfortunately, the technique has arguably produced 
more controversy than any other relation in astronony. 


The idea that that rotation rate of a spiral galaxy might be use- 
ful as a distance indicator goes back to Opik in 1922, but the 
first modern formulation was presented by Mort Roberts in the 
1960’s. Consider a system where some H I is distributed signifi- 
cantly outside the galaxy’s main mass concentration. According 
to Kepler’s law, the period of rotation of this H I gas is related 
to galactic mass by 
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where we assume a circular orbit and r is the distance of the gas 
from the galactic center. The rotational velocity of this gas is 
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or, since the radius is related to the observed angular size by the 
distance r = d0 
Miot « Od v2 (15.04) 


Now suppose we group galaxies together such that, for a class 
of galaxies, the ratio of hydrogen gas to total galactic mass is a 


constant. In other words, (g) = My 1/Mtot. Thus 
My 1 « Odv2 (15.05) 


Since the mass of H I present in a galaxy is related to the total 
flux emitted at 21 cm 


Muy 1 « d? fa «x Od? (15.06) 
where fo, is the observed 21 centimeter flux. Consequently, 
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Once could therefore observe nearby galaxies to determine the 
constant of proportionality, and then use the relation between 
total 21 centimeter flux and rotation velocity to derive galactic 
distances. 





Note here that radio telescopes like Arecibo do not have high 
spatial resolution. Therefore, unless the galaxy is nearby (i.e., 
has a large angular size) or is observed with a large array, it is 
difficult to determine the motion of the gas at particular locations 
in the galaxy. However, you can measure the total Doppler width 
of the H I emission, and the width of the emission line reflects 
the galaxy’s rotation. This is the way large samples of galaxies 
are investigated. 


(Roberts’ formulation for a luminosity-line-width relation was 
never seriously investigated, though it did prove to be a motivat- 
ing factor for today’s Tully-Fisher Relation.) 
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In most radio telescope observations, galaxies are unresolved (or 
just marginally resolved). 
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The Tully-Fisher Relation 


Perhaps the most important relation in extragalactic astronomy 
is the luminosity-line width relation for spiral galaxies, otherwise 
known as the Tully-Fisher relation. Its existence not only has 
implications for the extragalactic distance scale, but also for the 
structure of galaxies, the distribution of dark matter, the forma- 
tion of galaxies, and large-scale structure in the universe. 


Once again, let’s start with the idea that the rotation rate of a 
galaxy is somehow related to its mass. From above 


Miot x v2 (15.08) 


Now, let’s generalize this equation by saying that the quantity 
M /r is proportional to galactic luminosity, and that the coeff- 
cient of ve need not be 2. Thus 


Lau (15.09) 
or, in terms of magnitudes 
M =alogu. + b (15.10) 


Or, if we note explicitly that ve is not observed, but the Doppler 
line width (W) is, then 


M=alogW +b (15.11) 


This is the modern form of the Tully-Fisher relation. Once cal- 
ibrated (via nearby galaxies with known distances), a galaxy’s 
observed line-width can be used to determine its absolute mag- 
nitude. A comparison with apparent magnitude then yields the 
distance modulus. 


Although simple in theory, the details of the observations can 
be tricky. First, the observed rotation rate (or line width) of a 
galaxy depends on the galaxy’s orientation to the line-of-sight. 
Specifically, the observed rotation rate is related to the true ro- 
tation rate via the galaxy’s inclination, 


Cirus = Usps) Sid (15.12) 


so an estimate of inclination is a must. (Note that for highly 
inclined galaxies, sini ~ 1, so that the correction to line-width 
is small for edge-on systems.) In addition, even in the case of 
a perfectly face-on system, the emission-line widths have finite 
width, due to the z-motions in a galaxy. So one needs to know 
how important this component of “turbulence” is. 


In addition, galactic internal extinction will affect the measure- 
ment of the total apparent magnitude. The effect is largest at 
short wavelengths (say, B) and least in the infrared. But it is 
also inclination dependent. Highly inclined galaxies will have a 
larger fraction of their light extincted, since more of the galaxy’s 
luminosity must past through the dust on its way to the observer. 


The First Tully-Fisher Distances 


[Fisher & Tully 1977, Comm. Astrophys., 7, 85] 
[Tully & Fisher 1977, Astr. Ap., 54, 661] 
[Sandage & Tammann 1976, Ap. J., 210, 7] 


The first paper detailing the Tully-Fisher relation appeared in 
1977; it concluded that the distance modulus to the Virgo Cluster 
was (m—M) = 30.8 (14.4 Mpc). Sandage and Tammann quickly 
published a rebuttal paper, which appeared in 1976: using the 
same data, and applying the same technique, they concluded that 
Virgo’s distance modulus was (m— M) = 31.45 (19.5 Mpc). Why 
the difference? 
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Fig. 5 (a) Absolute magnitude — global profile width relation produced 
by overlaying Figure 3 on Figure 1, adjusting Figure 3 vertically to 
arrive at a best visual fit with a distance modulus of uo =30"6+0™2 
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Determining Inclination 


To first order, one can estimate the inclination of a spiral galaxy 
by assuming that the galaxy is a infinitely thin disk. When 
viewed face-on, the observed axis ratio, b/a = 1; when viewed 
edge-on, b/a = 0, so under this approximation, 


cosi = b/a (15.13) 


Unfortunately, galaxies are not infinitely thin, so (15.13) breaks 
down at high inclinations. A better estimate is obtained by mod- 
eling a spiral galaxy as an oblate spheroid, with axis ratios a:a:c, 
where c < a. It’s a relatively simply geometry problem to derive 
the formula 


cos? i = 


(b/a)? — a? 
1 — Q? 
where a = c/a is the intrinsic axis ratio, and b/a is the observed 

ratio. For Sc galaxies, œ is thought to be ~ 0.13. 


(15.14) 


The derivation of the above equation is left as an exercise to 
the interested student, but to get you started, consider that the 
equation of an ellipse can be written in parametric form as 





x = acos ô y = csin (15.15) 





To relate the size of the true minor axis to the size of the observed 
minor axis, consider that the latter is defined by the location 
along the ellipse where the slope is tangent to the line-of-sight. 
In other words, you need to know the value of 0 where 
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The length of the observed minor axis is then the vector from 
the center of the ellipse to the tangent point, projected onto the 














plane of the sky (which is defined as being perpendicular to i) 
b = (acos@, csin 0) - (1, = tan) (15.17) 


After some algebra, the result is equation (15.14). 


The above analysis, of course, requires that the observer know 
the true axis ratio, a. This can be estimated by observing a large 
number of galaxies, and seeing how the observed distribution of 
b/a agrees with the prediction. For an ensemble of randomly 
oriented galaxies, the inclinations should be distributed as sin 7 di. 
So 

(b/a)? — a’ 


1 — a? 


cos? i = 


(15.18) 


implies 


sinidi = (b/a) {[(b/a)? — a?] [1 - a?) }™ d(b/a) (15.19) 
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Fic. 7. Histogram of corrected minor to major diameter ratios for the gal- 
axies in our sample. The dashed line is the distribution function given by Eq. 
(19) with g=0.13. Curve and histogram are arbitrarily scaled so that the 
areas under them for b/a<0.7 are the same. 


[Giovanelli et al. 1994, A.J., 107, 2036] 


Inclination Errors 


Inclination measurements are critical for the success of a Tully- 
Fisher analysis, and any error in 7 will propagate directly into 
the derived distance. Consider a galaxy with observed axis ratio 
b/a. This ratio, naturally, will have a small amount of measure- 
ment error, Obja. Let’s see how this error propagates into the 
estimation of the rotational velocity. For simplicity, we’ll assume 
spiral galaxies are infinitely thin disks (the mathematics and re- 
sult for the oblate spheroid is extremely similar). If € = b/a, then 
cosi = €, and from (15.12), the true rotational velocity is 


Ut = Vobs/ Sin i (15.20) 


SO 


fe gs AE (15.21) 


(1—cos2i)/?7_ (1 — 2)” 


According to the rules of error propagation, 


GC Ho. a (15.22) 
Ut Vobs Oc 


Doing the math yields 


€ 
Ov, = OeUobs {Sa} (15.23) 


This is similar to the result one obtains for a non-zero value of 
Q, 


€ 
Ov, = Ogee = ery aaa} (15.24) 


The interesting property about this function is that the term in 
brackets is extremely dependent on inclination: at high inclina- 
tion, errors in inclination make very little difference, but at low 
inclination, a small uncertainty in b/a means a large error in the 
derived rotational velocity. 
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Errors in inclination will also affect the correction for galaxy 
luminosity (due to extinction). However, this uncertainty does 
not become terrible until very high 2 > 80° inclinations. 


Tully-Fisher and Internal Extinction 


[Holmberg 1958, Medd. Lunds Obs. Ser. 2, No. 136] 
[Tully & Fouque 1985, Ap.J. Supp., 58, 67] 


Dust can have an affect on the derived luminosities of galaxies. 
The zeroth order correction to this can be computed by assuming 
the dust layer is a simple slab parallel to the disk of the galaxy. 
In this case, the extinction is proportional to the path length 
through the dust, which, for an infinitely thin slab is 


A; = qseci = > A; — Aj=o = qseci — qsec(0) = q(sec i — 1) 
(15.25) 























Fits to Sc galaxies suggest, q ~ 0.28 mag (in the B-band). A 
better approximation is to say that the dust lies in a layer that 
cuts through the center of the galaxy. In other words, half the 
galaxy is extincted, and the other half (above the dust) is not, 
t.e., 

L=0.5L+0.5Le "8" (15.26) 


SO 
A; = —2.5 log {0.5(1 + e7" 5t) } 4527) 





Or, one can do even better: one can assume the dust sheet has a 
finite thickness. In this case, a fraction of the light f, is in front 
of the dust, a fraction f is behind the dust, and the remainder 
is intermingled. A decent fit is obtained using 


A; = —2.5 log ‘fa + e~"8ee") + (1 — 2f) (——)\ 


K Sect 
(15.28) 
with «x ~ 0.55, f ~ 0.25, and 7 is never given a value greater than 
80°. 


Alternatively, one can take an completely empirical approach. 
One can observe large samples of galaxies, divide them into 
redshift bins, and look for magnitude trends with size, surface 
brightness, and/or velocity line-width. Then you can just cor- 
rect the data using a value AM that depends on inclination (and 
bandpass). 


Using Tully-Fisher 


In practice, the calibration of the Tully-Fisher is fairly straight- 
forward. 


e Choose a cluster of galaxies that are presumably all at the same 
distance 


e Observe a complete sample of these objects at optical wave- 
lengths (to determine their relative magnitudes), and at radio 
wavelengths (to measure their H I Doppler line-width). In gen- 
eral, to minimize the effects of internal extinction and variations 
in the present star-formation rate, it is best to do the optical ob- 
servations in the red or infrared. (That way, a few bright blue O 
and B stars won’t affect the galaxy’s total luminosity that much.) 
Presently, CCD J-band observations are most common, but H- 
band aperture photometry has also been done extensively. ) 





e Plot magnitude versus the log of the line-width, and determine 
the slope of the Tully-Fisher relation. 


e Measure the magnitude and Doppler line-widths for a few 
galaxies with known distances (and therefore known absolute 
magnitudes). Use these absolute magnitudes to determine the 
zero-point (intercept) of Tully-Fisher relation. 


The Hubble Space Telescope was named that since one of its prin- 
ciple projects at inception was to measure the Hubble Constant 
to 10%. This task fell to the HST Distance Scale Key Project. 
Their measurement method was to determine the Cepheid dis- 
tances to a dozen or so local, high-inclination spiral galaxies and 
calibrate the zero-point of ‘Tully-Fisher. 





Malmquist Bias 
[Malquist 1922, Lunn Medd. Ser. I, 100, 1] 


An extremely important issue in all astronomical observations 
is that of biases. Even the most complete survey is subject to 
biases of one type or another that can affect the result. One of 
the most famous (or infamous) of these biases is Malmquist bias. 
This bias muddled the issue of Tully-Fisher distances for almost 
three decades. 


To understand what Malmquist bias is, let’s consider the original 
problem that Malmquist set out to investigate — the luminosity of 
F-stars in the Galaxy. Suppose we are dealing with a magnitude 
limited sample of stars, i.e., all objects brighter than a certain 
apparent magnitude are being analyzed. Let’s also suppose that 
not all F-stars are exactly alike, but instead that there is a Gaus- 
sian distribution of absolute magnitudes, with mean magnitude 
Mp and dispersion, ø. In other words, the true luminosity func- 
tion of stars is 





(M) = Poe (M—Mo)"/20° (16.01) 


Now let p(r) be the space density of stars, and A(m) be the 
observed luminosity function, i.e., 


AGS J “aU (16.02) 
A(m) = [ Boe (M—Mo)"/20" or) r°dr (16.03) 


Now for our analysis, we will need the derivative of this function 
with respect to apparent magnitude. Since, from the definition 
of distance modulus, 


— = — (16.04) 


we can write 
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Now, if we divide by A(m), we get 





dA Mo 1 f, M®(M 2d 
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dm g- 2 fo ®(M)p(r) r2dr 

But note that the last term is, by definition, the mean absolute 

magnitude of the sample. Thus 


Mo — (M) = 0? = (16.07) 
In other words, the mean absolute magnitude of the stars in the 
sample does not equal the intrinsic mean for the set of stars. 
In particular, the difference is proportional to the dispersion of 
the sample, squared. (If the intrinsic dispersion in the sample is 
large, you can be off by quite a bit!) 


Malmquist bias comes about because of the artificial cutoff im- 
posed by a magnitude limited sample. If the apparent magnitude 
cutoff is m, then more-or-less all objects out to a distance 


gate (16.08) 


are being observed. However, some exceptionally bright objects 
at larger distances will make it into the sample, while some ex- 
ceptionally fainter objects that are closer than d will not. Thus, 
you are biased towards brighter objects. 


Note that this is not just a minor inconvenience. Suppose you 
wanted to derive the distance to a galaxy by observing its F-stars. 
(They’d be too faint, but never mind that.) If you did not include 
the effects of Malmquist bias, then you would overestimate the 
true absolute magnitudes of F-stars in the Milky Way, and, since 
(m—M) = 5log(d/10), you would underestimate the distance to 
the galaxy. (And, since the Hubble Constant is V/D, you would 
therefore overestimate Ho.) 





Malmquist Bias — An Example 


Let’s look at a specific example of Malmquist bias, where the 
assumed space density of objects is constant. The total number 
of objects within a distance, d, will then be 


d 
1 
N(d) = J pr°dr = ged (16.09) 


For stars of absolute magnitude M, d is related to apparent mag- 
nitude, m, via (16.08), so 
1 1 
N = = 100 6(m-M)+3 a 1.382(m—M)+6.908 16.10 
(m) = 5p “pe (16.10) 
In other words 
N(m) x gsm (16.11) 


But note that N(m), the total number of objects brighter than 
m, is just the integral of the observed luminosity function down 
to apparent magnitude m, 


NGA= J Aada (16.12) 
SO _ 
A(m) = aN) oc 1.382607" (16.13) 
i Ce A _ 1 389 (16.14) 
dm 


This means that the objects in the sample will be brighter than 
the norm by a factor that goes as 


(M) = Mo — 1.38207 (16.15) 


So, in order to derive the correct absolute magnitude, you need 
to know the dispersion of the distribution and have some idea of 
the true space distribution of objects. 


Malmquist Bias and the Determination of c 


The situation is even worse than it sounds. Let’s use the same 
magnitude limited sample of objects, and try to estimate o. To 
do this, we start with (16.05) 
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and compute the second derivative 
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Collecting terms gives 
,d7 n A(m 
ie nan) X = M —2Mo(M) + (M?) — o?” (16.18) 


or 


g fı + | = (Mo — (M)? (16.19) 


But note: the right hand side of this equation is just the ob- 
served dispersion in magnitudes! It is not the same as the true 
dispersion. 
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Fig.l. Upper panel: Monte Carlo distri- 
bution in distance and absolute magnitude 
of 500 galaxies within 38 Mpc. Constant 
space density and a mean absolute magnitude 
of <M>=-18™ with a Gauss standard deviation 
of oy=2™ are assumed. Lower pannel: The 
same sample cut by an apparent-magnitude 
limit of m=13™0. Note the increase of the 
galaxian luminosities with increasing 
distance and the small effective (observable) 
scatter oj), within individual distance 
intervals. 


Malmquist Bias and the Tully-Fisher Relation 
[Teerikorpi 1997, A.R.A.A., 35, 101] 


The previous derivations for Malmquist bias apply if the intrinsic 
luminosity function of objects is a Gaussian with some dispersion. 
In the case of the Tully-Fisher relation, we are not working with 
a single mean magnitude, but with a relation between magnitude 
and rotation velocity. Under these circumstances, the mathemat- 
ics of correcting for Malmquist bias is more complicated, but, in 
general, the result is the same. The effect of Malmquist bias is to 
change the slope of the Tully-Fisher relation by an amount that 
goes as the square of the dispersion about the mean line; the more 
the scatter, the steeper the apparent slope and the more you un- 
derestimate distance. As a result, Tully-Fisher arguments often 
focus on the amount of scatter in relation. For instance, one ob- 
serves a set of galaxies in a cluster, and measures the slope and 
scatter in the Tully-Fisher relation. Is the scatter intrinsic to 
the Tully-Fisher relation, or is there depth and/or foreground- 
background contamination in the cluster? Depending on your 
point of view, you can use the observations to support a small or 
large Hubble Constant. 





(The paragraph above may be debated by some. In particular, 
one argument goes that if you take a velocity-line width limited 
sample of galaxies, or select objects by H I line width (instead 
of optical magnitude), you can avoid the issue of Malmquist bias 
completely. But these formulations have their own statistical 
problems. ) 


In the 1990’s, the debate about the Hubble Constant largely 
came down to an argument about the intrinsic scatter in the 
Tully-Fisher relation. In general, high Hubble Constant groups 
claimed that the scatter about the Tully-Fisher line was about ~ 


0.3 mag; low Hubble-Constant people said the scatter is closer to 
~ 0.7 mag. Since distance goes as the square of ø, this translates 
into a ~ 0.6 mag difference in distance modulus, or almost a 
factor of 2! 


The answer seems to be that the Tully-Fisher relation does in- 
deed have a small amount of cosmic scatter. Observations in 
the Ursa Major Group show a dispersion of about the mean of 
o ~ 0.35 mag. The larger (o ~ 0.8 mag) dispersion in the Virgo 
Cluster seems to be due to the intrinsic depth of the cluster: 
Virgo is elongated along our line-of-sight. Note that the disper- 
sion in the ‘Tully-Fisher relation is small, and that the rotation 
rate of spirals is largely defined by the mass in their dark matter 
halos, the precision of the relation says a lot about the uniformity 
of the formation of disk systems. 


Note also that the scatter in the relation depends on color: the 
B-band relation has significantly more scatter than the infrared 
law. This is most likely due to the effect of dust, and the deem- 
phasizing of blue stars. (Observations in the blue record the 
young stars, but IR measurements probe the bulk of the stellar 
mass. ) 
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Fic. 11—B-, R-, I-, and H-band Tully—Fisher relations for the Local Calibrators (top), Ursa Major cluster members (middle), and Virgo cluster 
members (bottom). It is apparent from the figures that the slope of the relations increases going to longer wavelengths and the dispersion decreases. The 
variation in slope is thought to arise from the differing contributions to the observed bandpass made by greater fraction of young stars found in the 
lower-luminosity systems. The smaller dispersion at longer wavelengths is likely due to a reduction in the sensitivity to these effects, as well as those 
expected from extinction variations. Note the much larger dispersion found for the Virgo cluster data. 


